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Using the degenerate double exchange Hamiltonian with on-site Coulomb interactions we show 
that there is an instability towards the displacement of the manganese ion in Lai_a;Cai,Mn03. The 
result is a dipole moment due to the charge disproportionation and the lattice distortions that gives 
rise to the predicted magnetic ferroelectric phase for dopings in between x = A and x = 5 Q. 
The instability is stabilized by the phonons of the lattice, resulting in a stable configuration of the 
degrees of freedom of the system. 



Recently, manganites Ri-ajA^jMnOs (R=:rare-earth 
ion, A=divalent substituent), with their particular colos- 
sal magnetoresistance, have been the subject of intense 
study E S II & II IE 13 H- The internal degrees of 
freedom, charge, orbital state, and spin, vary with the 
composition x, resulting in different properties. Inter- 
estingly, these compounds very rarely are both magnetic 
and ferroelectric, which would be very desirable for tech- 
nological applications. 

It is well known that ferroelectric perovskite com- 
pounds originate in the displacement of the mostly non- 
magnetic transition metal ion from the center of the oc- 
tahedron The result is a stronger covalent bond be- 
tween the transition metal ion and one of the surrounding 
oxygens. The existence of a ferroelectric state in under- 
doped manganites due to the charge order that results in 
a lack of inversion symmetry has also been shown J| . We 
show that the lattice distortions due to the displacement 
of the magnetic manganese ion, in addition to the previ- 
ously discussed interaction between charge, orbital, and 
spin ordering in underdoped manganites [l|, [lOj , further 
lower the energy of the ferroelectric state, and give rise 
to an electric dipole moment. The result is the desired 
ferroelectric magnet compound. 

The main complication arises from the anisotropy of 
the orbital states of the manganese ion. The fivefold de- 
generacy of its d orbital states is partially lifted by the 
crystal field into twofold degenerate orbital states, d^^-y-^ 
and d322_r2, called Cg, and threefold degenerate orbital 
states, dxy, dyz, and dzx, called t2g. The degeneracy of 
the Cg states is further lifted by the lattice distortions, or 
by the on-site Coulomb interactions. The former lift the 
degeneracy completely by lowering the symmetry of the 
cubic crystal which lowers the energy of the d^z^_r2 state 
while raising that of dj.2_j^2. The latter make the system 
quasidegenerate because the degeneracy is only lifted af- 
ter either of the two orbital states is first occupied. After 
choosing an orbital basis we maycharacterize the orbital 
states of every ion by an angle 

The electronic configuration of each ion in the array 
may also vary. The Mn^+ ion in the parent compound 
of interest, LaMnOa, has an open shell configuration in 
which three electrons are in the t2g orbitals and one in ci- 
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FIG. 1: The lattice with the corresponding unit cell for the 
intermediate phase (j> — 7r/4. The CE phase corresponds to 
4> = 0, and the Zener polaron to = 7r/2. The Zener polaron 
phase is made up of dimers, with blue and red dimers AFM 
coupled within a diagonal. The CE phase is composed of 
AFM coupled zig-zag chains. The intermediate phase shows 
the superposition of the Zener polaron and the CE phase. 



ther of the two eg orbitals. The t2g electrons are localized 
due to a stronger stabilization of the crystal field with re- 
spect to the Cg orbitals and a lower hybridization with 
the oxygens 2p states ^3 ■ When doping with a divalent 
ion that introduces holes, for instance Lai_2;Caa;Mn03, 
some sites change to Mn''+ with a closed shell configura- 
tion in which the Cg orbitals are empty. 

There are three interactions that reproduce the prop- 
erties of the solids in discussion: the double-exchange, 
the superexchange, and the on-site Coulomb interaction. 
The double-exchange, associated with the hopping am- 
plitude t, leads to the formation of bands and to the con- 
duction of electrons. The superexchange, with the energy 
scale J, is a magnetic interaction between localized spins 
of neighboring sites. The on-site Coulomb interaction is 
between electrons in different Cg orbital states in the same 
ion. Its energy scale is U, and in the solids of interest we 
find U >> t, which tends to localize the electrons. 

We consider two limiting cases of Lai-jjCa^^MnOa: the 
Zener polaron phase at a; = 0.4 doping, and the CE phase 
at X = 0.5. Both phases are known to be insulating. As 
shown in Fig. ^ the CE phase is made of zig-zag chains 
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with ferromagnetic spin order within the chain, and that 
are antiferromagnetically coupled with the neighboring 
chains. This phase has the checkerboard charge order 
in which sites alternate between Mn'^+ and Mn''+. The 
unit cell has four sites, with one Mn ion each. The orbital 
order is d^^2_j.2 and dx2_y2 in the corner sites (CI and 
C2), d3a;2_^2 in the horizontal bridge (B2), and d3j,2_^2 
in the vertical bridge (Bl). The Zener polaron phase 
is made up of dimers with equivalent sites (B1=C1 and 
B2=C2). The spins within a dimer are parallel, while 
spins in between neighboring dimers with the same crys- 
tallographic direction couple antiferromagnetically. Per- 
pendicular dimers have perpendicular spin orientations 
with respect to each other. The charge in a dimer is de- 
localized in between the equivalent sites, which results in 
a bond-centered charge order. 

The Hamiltonian of the system is given by 

<ij> .(T(t' 

+ Jaf Y^'^-^i + U J2 ' 



nlnl , (1) 



where j,j are site labels, and tr, cr' denote Cg orbital 
states. The sum < i] > is over nearest neighbor sites. 
The on-site Coulomb interaction contains the density per 
eg orbital per site, nf, and Jaf is the antiferromagnetic 
superexchange interaction between (classical) localized 
iig spins, Si. The kinetic term with the hopping integral 
t contains the fermionic creation and annihilation oper- 
ators, cj^ and Cj^o-', and the symmetry matrix elements 
Ff/ that take into account the overlap of the orbitals 
involved and the crystallographic direction of the hop- 
ping. We take into account the overlap of the spin wave 
functions 0|, (xi|xj), which depends on the relative di- 
rection between the two spins, | {xi\xj) I = cos((/)y72). 
For parallel spins the hopping is largest, and antiparallel 
spins cannot hop. 

We consider the zero temperature limit with strong 
ferromagnetic Hund's coupling between localized t2g and 



itinerant eg electrons, J/z/t 



The antiferromag- 



netic interaction becomes J AF^^ij--, jS'p cos(/)y, where 



the angle denotes the relative orientation between 
spins. The localized spins have jS'p = 9/4. With re- 
spect to the symmetry matrix, we choose our orbital ba- 



sis states to be a 
case we have 
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that takes into account the overlaps of the Cg orbitals in- 
volved when hopping along the y-axis . The hopping 
along the x-axis has a phase with respect to the y-axis 
which results in = ^%>//^, = "VS. After 

choosing a basis we may refer the orbital states of every 



ion, ai and Pi, to the reference basis with the transfor- 
mation 
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cos(|)|q; > -|-sin(|)|/3 



- sin(|)|a > -I- cos(| 
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Finally, we take the on-site Coulomb interaction to be 
U/t = 4 and use the mean-field approximation of the 
densities < nf >mf—< c\^^Ci,a >. 

After decoupling the on-site Coulomb interaction, 
nfnf < > nf + < < > - < < >< nf >, 
we first calculate the charge order (CO), orbital order 
(00), and magnetic order (SO) for different fixed val- 
ues of the doping varying the spin canting angle from 
(pij = (CE phase) to (j)ij = 7r/2 (Zener polarons). Since 
we consider the relative orientation between spins, (jjij, 
to change equally, in what follows we use a general spin 
canting angle, (j). We also consider that the relative angle 
does not change for some neighboring spins, as seen in 
Fig. n The anisotropy of the system is taken into ac- 
count by rotating the orbital states per site to our work- 
ing bases, a and /?, with the transformation 



,/3 



<• sin2(e,/2)-f nf' cos2(e,/2) 
n^* cos2(e./2) -fnf- sin2(e,/2) 



(4) 



before solving numerically. The self consistency equa- 
tions, 

= tan-i(2 <c^c.,^>/«-nf)) (5) 
nf' = <cos2(C,/2)-H<4„c,/3>sin(eO 

+ nf sin2(e,/2) 
nf = nf sin2(^,/2)- < cl^c,j3 > sm{^^) 

+ nf cos2(e./2), 

enforce that the rotation of the calculated densities, nf, 
from the working basis back to the on-site basis results 
in a diagonal on-site Coulomb interaction. 

Now comes the totally new result: we also take into 
account the displacement of the manganese ions by con- 
sidering small changes in the hopping amplitude, t — > 
t + c5, where 5 is the distortion of the lattice and the con- 
stant c = 3.5 is obtained from Ref. We first fix the 
spin canting while varying 5. Regarding the orbital states 
of the CE phase, shown in Fig. [21 the effect of the dis- 
tortions is to rotate the orbital states at the corner sites 
towards the bridge site orbital states making them more 
equivalent. For the equivalent orbital states of the Zener 
polaron phase there is no change, which shows that the 
orbitals are independent of the distortion. Considering 
the energy up to linear order in S and without the elastic 
energy quadratic in S we find an almost constant energy 
for the CE phase, while for the Zener polaron phase the 
energy decreases linearly. Hence, there is an instability 
towards the distortion of the lattice. 
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FIG. 2: Orbital order for fixed doping (top row), and energy 
without elastic energy contribution (bottom). The left col- 
umn corresponds to the CE phase, and the right to the Zener 
polaron phase. 



03 

^5 



0.04 



0.02 




FIG. 4: Charge disproportionation (green), and lattice distor- 
tion (red) as a function of doping. The Zener polaron phase 
has no dipole moment due to Aq = 0, and the CE phase since 
5 = 0. 




FIG. 3: The nonzero value of Fo{x) indicates an instability to- 
wards lattice distortions. Only the CE phase does not present 
the instability. 



In order to find the instability, we calculate the deriva- 
tives of the energy with respect to the distortions for the 
expansion 



E — En 



dE ^ 



(6) 



The rate of change of the energy with respect to the dis- 
tortion relates to the rise of ferroelectricity. Introducing 
5 ~ S/a, the energy we find has the form 



E{d,x) = (^FoixyS + Eo{x)y, 



(7) 



where Fq{x) is a force, and Eq{x) is the ground state en- 
ergy without distortions. The force, F{S,x), is obtained 
by taking the derivative of the energy with respect to 
the distortion. We find that Fo{x) is the instability of 
the undistorted lattice towards lattice distortions. The 
energies without the elastic energy contribution decrease 
linearly as a function of the distortions. Thus, we con- 
clude that the addition of the elastic energy will stabilize 
the system. As seen in Fig. 01 the CE phase has no in- 
stability, while for the Zener polaron phase we see a large 
instability. 

Subsequently, we add the elastic energy term, 
{l/16)Dd^, where D is the elastic constant, in order to 



find the minimum energy Eo{S, x). Defining D — Da/lQt 
the energy now has the form 



E{5,x) = {bp -f Fq{x)~5 -f Eq{x)^ t. 



(8) 



Its minimum is given by the value of the distortion, for 
which F{5{)^x) — 0. Taking the derivative of the energy 
with respect to the distortion we find 



(5o 



Fq{x) 
2D 



(9) 



The elastic constant D can be calculated for the Zener 
polaron phase with D — — _Fo(2;)/2(5n, with 6o « 0.042 
obtained from experimental results fSj, where a is the 
undistorted lattice spacing between sites. Hence, we find 
a value of D = 16.62, which we consider to be indepen- 
dent of the doping. Using t = 0.622 eV from Ref. [l^ and 
a = 3.947 A from Ref. ^ we recover D = 10.61_eV/A2, 
which is in agreement with other calculations |17l |. 
Therefore, the dipole moment, 



p = SAq, 



(10) 



is zero for the Zener polaron phase and the CE phase. 
Nonetheless, the cause is different: in the case of the 
CE phase, it is due to the lack of distortions, and for the 
Zener polaron phase due to a vanishing charge dispropor- 
tionation, as shown in Fig. ^ Strikingly, the in between 
phases have a nonzero electric dipole moment, shown 
in Fig. that is directly proportional to the on-site 
Coulomb interaction: Aq increases as U increases. This 
counter intuitive result states that the on-site Coulomb 
interaction promotes a metallic state instead of localizing 
the electrons for dopings in between 0.4 and 0.5. 

Finally, we extend our calculation to also include vari- 
ations of the spin canting. The result is shown in Fig. |H| 
for the Zener polaron phase and CE phase. For the CE 
phase we find that the spin canting is independent of the 
orbital and lattice distortion degrees of freedom. The 
Zener polaron phase shows only a small range where the 
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FIG. 5: Electric dipole moment p as a function of doping 
X. The maximum value of the dipole moment is p,nax ~ 7.5 
meA at around x ~ 0.44 
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FIG. 6: Orbital order (top), spin canting (middle), and en- 
ergy without elastic energy contribution (bottom). The left 
column is for the limiting case x = 0.4 and the right for 
X = 0.5. For the range in which the energy changes linearly 
with the distortion the spin and orbital degrees of freedom do 
not change. 

spin canting and the orbital order are independent of 
the lattice distortions. Using again the lattice distor- 
tions 5 « 0.042 for the Zener polaron phase, we see two 
effects: firstly, the energy without the elastic energy de- 
creases linearly with the distortions up to this value of 
5, and secondly, the orbital order and spin canting now 
have small corrections. Interestingly, the effect on the or- 
bital order is to make the two sites within a dimer slightly 
unequivalent. For larger distortions that go beyond the 
stable point, the energy changes its dependence on the 
distortion, and we find a larger rate of change of the spin 
canting and orbital states of the underdoped states to- 
wards a CE phase configuration. The elastic constant 
becomes D — 14.21, which means that the energy of the 
phonons decreases with the spin canting. In this case we 
recover D = 9.07 eV/A^. For other dopings within the 
range here discussed we expect the same behavior. 



We conclude that in underdoped manganites, such as 
Lai_2;Caa;Mn03, there is an instability towards lattice 
distortions for the states in between 0.4 and 0.5 dop- 
ing. This gives rise to a dipole moment that results in 
a ferroelectric state. Counterintuitively, stronger on-site 
Coulomb interactions result in a larger dipole moment, 
promoting ferroelectricity. That charge order can exist in 
a ferromagnetic region has already been experimentally 
discussed We have shown that the spin and orbit 
degrees of freedom are coupled to the distortion of the 
lattice. It has been proposed that including the oxygens 
in the calculations results in the CE phase being made 
up of oxygen stripes 4] . Interestingly, in our findings the 
magnetic correlation between dimers in the same stripe 
is always FM and does not change for the dopings be- 
ing discussed. We raise the question of whether these 
two results are related. We do not find any agreement 
with the results in which a different configuration of the 
CE phase was proposed 0. In addition, our results do 
not agree with the claim that Zener polarons are formed 
within the zig-zag chains of the CE phase and that they 
are important for understanding the localization of elec- 
trons in half-doped manganites Future work should 
include nonzero temperature calculations. 
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